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We calculate the thermodynamic quantities (correlation functions (SoS n ), correlation length £, 
spin susceptibility x, and specific heat CV) of the frustrated one-dimensional spin-half J1-J2 Heisen- 
berg ferromagnet, i.e. for J2 < 0.25|Ji|, using a rotation-invariant Green's-function formalism and 
full diagonalization of finite lattices. We find that the critical indices are not changed by J2, i.e., 
X = yoT~ 2 and £ = xqT^ 1 at T — > 0. However, the coefficients j/o and a;o linearly decrease with 
increasing J2 according to the relations yo — (1 — 4J2/|Ji|)/24 and Xo = (1 — 4J2/|Ji|)/4, i.e., 
both coefficients vanish at J2 = 0.25|Ji| indicating the zero-temperature phase transition that is 
accompanied by a change in the low-temperature behavior of \ (£) from \ x T~ 2 (c; oc T -1 ) at 
J 2 < 0.25|Ji| to x oc T" 3/2 (£ oc T" 1/2 ) at J 2 = 0.25|Ji|. In addition, we detect the existence of 
an additional low-temperature maximum in the specific heat when approaching the critical point at 
J 2 =0.25|Ji|. 

PACS numbers: 



I. 



INTRODUCTION 



Low-dimensional quantum magnets represent an ideal 
playground to study systems with strong quantum and 
thermal fluctuations X In particular, much attention has 
been paid to the one-dimensional (ID) J1-J2 quantum 
Heisenberg model, which may serve as a canonical model 
to study frustration effects in low-dimensional quantum 
magnets. Although this model has been studied fre- 
quently (see Ref. and references therein), the model 
deserves further attention to detect unknown features of 
this quantum many-body system, especially in the case 
of ferromagnetic nearest-neighbor (NN) interaction Ji < 
()4— — — — — — From the experimental side, recent stud- 
ies have demonstrated that edge-shared chain cuprates 
represent a family of quantum magnets for which the 
ID J\ — J2 Heisenberg model is an appropriate start- 
ing point for a theoretical description. Among others, 
we mention L1VCUO4, LiCu2C>2, NaCu2C>2, L^ZrCuO,}, 
and LizCuOzr^^^^^^^^^^- which were identi- 
fied as quasi-lD frustrated spin-1/2 magnets with a ferro- 
magnetic NN in-chain coupling J\ < and an antiferro- 
magnetic next-nearest-neighbor (NNN) in-chain coupling 
J2 > 0. The Hamiltonian of their ID subsystems consid- 
ered in this paper is then given by 



H — J\ 2_^ SiSj + J2 2_^ SiSj, 



(1) 



(id) 



[U] 



where (i,j) runs over the NN and [i,j] over the NNN 
bonds. For the model (TJ]) the ferromagnetic ground state 
(GS) gives way for a singlet GS with spiral correlations 
at the critical point J 2 = 0.25|Ji|i2i22i 

The edge-shared chain cuprates have attracted much 



attention due to the observation of incommensurate 
spiral spin ordering at low temperature. Hence, in 
these compounds the antiferromagnetic NNN exchange 
J2 is strong enough to destroy the ferromagnetic GS 
favored by the ferromagnetic 3\. On the other hand, 
several materials that considered as model systems 
for ID spin-1/2 ferromagncts, such as Tetramcthylam- 
monium Copper Chloride ( TMCuC [(CH 3 )4NCuCl 3 ]) 
(Ref. [231 ) and p-nitrophenyl nitronyl nitroxide (p-NPNN) 
(Ci3Hi6N304)f2i might have also a weak frustrating 
NNN exchange interaction J2 < — 0.25Ji. Moreover, re- 
cent investigations suggest that Li2Cu02 is a quasi-lD 
spin-1/2 system with a dominant ferromagnetic J\ and 
weak frustrating antiferromagnetic J2 ~ 0.2 1 Ji|^ 

Although for J2 < — 0.25 Ji the GS remains ferromag- 
netic, the frustrating J2 may influence the thermodynam- 
ics substantially, in particular near the zero-temperature 
critical point at J2 = 0.25| J\\, The investigation of this 
issue is the aim of this paper. The study of the ID J1-J2 
Heisenberg model is faced with the problem that, due to 
the J2 term, neither the Bethe-ansatz solution nor the 
quantum Monte Carlo method is applicable. Hence we 
use (i) the full exact diagonalization (ED) of finite sys- 
tems of up to N = 22 lattice sites, and (ii) the second- 
order Green's-function technique^ that has been applied 
recently successfully to low-dimensiona l q uantum spin 



systems 



26.27.2N. 2!) 



For example, in Ref. 1271 . by compar- 



ison with Bethe-ansatz data it has been demonstrated 
that this method leads to qualitatively correct results for 
the thermodynamics of the ID Heisenberg ferromagnet in 
a magnetic field. As the most prominent feature, a field- 
induced extra low-temperature maximum in the specific 
heat has been found 2 ^ and characterized as a peculiar 



quantum cffcct i 27 ' 29 



FULL DIAGONALIZATION OF FINITE 
LATTICES 



II. 



Using Schulcnburg's spinpack (Rcf. |30T ) and exploit- 
ing the lattice symmetries and the fact that S z — J2i &f 
commutes with H, we are able to calculate the exact ther- 
modynamics for periodic chains of up to N = 22 spins. 
The comparison of results for N = 12, 14, 16, 18, 20, and 
22 allows to estimate the finite-size effects. The largest 
matrix which has to be diagonalized for N — 22 has 
29414 x 29414 matrix elements. 
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III. SPIN-ROTATION-INVARIANT 
GREEN'S-FUNCTION THEORY 



FIG. 1: NN (solid) and NNN (dashed) spin-correlation func- 
tion for J2 = 0,0.1, and 0.2, from top to bottom, calculated 
by RGM (lines) and ED (open symbols; N=20). 



To calculate the spin correlation functions and the 
thermodynamic quantities, wc determine the trans- 
verse spin susceptibility Xq ( w ) = ^{{S q ;SZ q ))u 
(here, ((...;... )) u denotes the two-time commuta- 
tor Green's functional) by the spin-rotation-invariant 
Green's- function method (RGM) . 25 ' 26 Using the equa- 
tions of motion up to the second step and suppos- 
ing rotational symmetry, i.e., (Sf) = 0, wc obtain 
u> 2 ((S+;Sl q )) u = M q + ((-§+; SZ q )} u with M q = 
([[S+,H] ,SZ q }) and -5+ = [[S+,H] ,H]. For the 
model ([1]) the moment M q is given by the exact expres- 
sion 



M q = -4 Y^ J nC n (1 - cosnq) , 
n=l,2 



(2) 



where C n = (S+S~) = 2(S§S*). The second deriva- 
tive — S£ is approximated as indicated in Refs. |25||261 
l^HilM That is, in -§f we adopt the decoupling 
5,j + Sj~SZ = a(Sj~ SZ)Sf~ + a(S^SZ)Sj~ , where in the 
case J 2 < — 0.25Ji with a ferromagnetic GS the ver- 
tex parameter a can be assumed in a good approxima- 
tion to be independent of the range of the associated 
spin correlators (see the discussion below). We obtain 
2 <?+ and 



-bg 



with 



<=Y. 



<4s+ 



xt M 



-\{S q ; S_ q )) 



Mr, 



«?-w 2! 



(3) 



J n Jm (1 - cosnq) [K ntm + 4aC„ (1 - cos mq)] 



n,m(=l,2) 



(4) 
where K n ,„ = 1 + 2a (C 2n - 3C„), K x , 2 = 2a (C 3 - Ci), 
and K 2 ,i = K\% + 4a (C\ — C 2 ). From the Green's func- 
tion {3J the correlation functions C„ = -^ J2 C q e lqn 
of arbitrary range n arc determined by the spectral 



theorem^ 



C q = (S+SZ q ) = ^[l + 2n(c q )} 



2uj„ 



(5) 



where n(u) q ) = (e w "/ T — l) 
the operator identity S^~S~ 



is the Bose function. By 
= h + Sf we get the sum 



rule Co = jf J2 q C q 
lim 



ccptibility x 
Xq M = \x\ 



- 4. The uniform static spin sus- 

q-^oXq, where \q = Xq (^ = 0) and 
(u>), is given by 






71=1,2 



n,??i(— 1,2 



TL Jri'Jri 



K n 



(6) 

The correlation length £ may be calculated from the ex- 
pansion of the static spin susceptibility around q = (see, 
e.g., Refs. [H and[2g|) x g = xl (l + £V) • The ferromag- 
netic long-range order, occurring in the ID model at T 
only, is described by the condensation term C (Ref . [2f 
according to C„ (0) = £ £^0) (M q /2u q ) e * 9 " + C- 
Equating this expression for n ^ to the exact result 



C„ #0 (0) 



the ratio M q /2uj q must 



(So<Wo) (°) = 
be independent of q, because C„^o is independent of n. 
This requires the equations K n ^ m (0) = [cf. Eqs. ([2| and 
Q], which yield a(0) = |. Then, wc get uj q (0) = f M g (0) 
and C = g, where the sum rule Co = 5 is fulfilled. In 
Eq. © we have A(0) = 0, so that x diverges as T — > 
indicating the ferromagnetic phase transition. 

Let us discuss the used assumption that the vertex 
parameter a is independent of the distance I. For that 
we consider an extended decoupling with four different 
parameters a/ (I = 1, ...,4) attached to the four corre- 
lators Ci appearing in w 2 [cf. Eq. ([J}]. At T = 0, the 
four equations if„ ]m = (n, m = 1,2) yield the solutions 
a;(0) = |. On the other hand, in the high-temperature 
limit all vertex parameters approach unity^ Because 
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FIG. 2: Uniform static spin susceptibility calculated by RGM 
(solid lines) and ED (open symbols; N=20) for J 2 = 0,0.125, 
and 0.2, from right to left, and by TMRG (filled symbols) for 
J-i = and 0.125 (Ref. |g). The inset shows the correlation 
length obtained by RGM for J 2 = 0, 0.125, and 0.2, from right 
to left. 



we have identical vertex parameters at T = and for 
T — > oo, we put ai = a in the whole temperature region, 
as was done above. 

To evaluate the thermodynamic properties, the corre- 
lators Ci (I = 1, . . . ,4) and the vertex parameter a have 
to be determined as numerical solutions of a coupled sys- 
tem of five non- linear algebraic self-consistency equations 
for Ci including the sum rule Cq ~ ^ according to Eq. |(5J). 
Tracing the RGM solution to very low temperature, we 
find that it becomes less trustworthy for J 2 approaching 
J 2 = 0.25|Ji|. Therefore, below we will present RGM 
results for J 2 < 0.2 1 J\\ only. 



IV. RESULTS 

Hereafter, we put |Ji| = 1. First we consider the NN 
and NNN correlation functions shown in Fig. Q] The 
RGM results agree qualitatively well with the ED data. 
Note that the difference between ED and RGM results at 
low temperature might be partially attributed to finite- 
size effects in the ED data. For larger temperature T > 1, 
the agreement becomes perfect. With increasing frustra- 
tion the correlation functions decrease, where the NNN 
and further-distant correlators decay much stronger than 
the NN correlator (interestingly, for J 2 = 0.2 the NNN 
correlator changes the sign at T w 1). This frustration 
effect is reflected in the correlation length £ depicted in 
the inset of Fig. O At T = 0, £ and the uniform static 
spin susceptibility x diverge due to the ferromagnetic GS. 
With growing temperature the decay of £ increases with 
increasing J 2 . As shown in Fig. [2J our ED data for \ are 
in excellent agreement with the results of the transfer- 
matrix renormalization-group (TMRG) study of Ref. (J 
and agree well with the RGM results. The susceptibility 







FIG. 3: xT 2 versus \/T calculated by RGM (solid lines) and 
ED (N = 20; dashed lines) for J 2 = 0,0.125, and 0.2, from 
top to bottom. For comparison we present also Bethe-ansatz 
data (open squares) for J2 = (Ref. |32T ) and TMRG data 
(filled triangles) for J2 = 0.125 (Ref. 0). The upper inset 
shows the coefficient yo — liniT^o XT 2 obtained by the RGM 
(filled squares) and ED (open circles) in dependence on J2 as 
well as a linear fit of the RGM data points (solid line) . The 
lower inset shows the coefficient y\ [cf. Eq. (J?}] obtained by 
the RGM (filled squares) and ED (open circles) in dependence 
on J2 as well as a quadratic fit of the data points (solid line). 



decreases with increasing J 2 , because this antiferromag- 
netic interaction counteracts the spin orientation along a 
uniform magnetic field. 

Next we investigate the critical behavior of x an d £ 
for T — > in more detail. To study the influence of the 
frustration on the critical behavior we follow Refs. [32| 
and |33| . The critical indices 7 for x an d v for £ can be 
obtained by analyzing the RGM data for — ^IcTTtT an< ^ 

- aloejr) for T ~* °- Wc nnd that 7 = 2 and 2/ = 1 are 
independent of J 2 for J 2 < 0.25. Going beyond the lead- 
ing order in T we know from Bethe-ansatz dat a 32 i 33 and 
from the renormalization-group technique 3 -^ that the low- 
temperature behavior of the susceptibility and the corre- 
lation length of the unfrustrated ID spin-1/2 Heisenberg 
ferromagnet is given by 

X T 2 = y Q + yiVT + y 2 T + 0(T 3 / 2 ) (7) 



and 



£T = x + x 1 Vt + x 2 T + O{T 3/2 ) . (8) 



Here we adopt this expansion suggested by the existence 
of the ferromagnetic critical point at T = 0, but with J 2 - 
dependent coefficients for the frustrated model (p}. To 
determine the coefficients yo and xo, in Figs. [3] and [4] we 
show the quantities %T 2 and £T versus y/T. Again wc 
find a good agreement of the ED for %T 2 with Bethe- 
ansatz and TMRG data down to quite low temperature. 
The RGM results for xT 2 and £T deviate slightly from 
the Bethe-ansatz and TMRG data for finite temperature. 
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FIG. 4: £T versus VT by the RGM (solid lines) for J 2 = 
0,0.125, and 0.2, from top to bottom. For comparison we 
present also Bethe-ansatz data (open squares) for Ji = 
(Ref. |33). The left inset shows the coefficient xo = liniT^o £.T 
obtained by the RGM (filled squares) in dependence on J2 as 
well as a linear fit of the RGM data points (solid line) . The 
right inset shows the coefficient x\ [cf. Eq. (|5J] obtained by 
the RGM (filled squares) in dependence on J2 as well as a 
quadratic fit of the data points (solid line). 
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FIG. 5: Specific heat obtained by RGM (solid lines), ED 
(open symbols; N=20) and TMRG (filled symbols; Ref. H) 
for J2 = and 0.125, from top to bottom. 



The behavior of the leading coefficients yo and xq and 
the next-order coefficients y\ and X\ can be extracted 
from the data for yT 2 and £T by fitting these data to 
Eqs. (O and ([5]). For the RGM we use data points up 
to a cut-off temperature T = T cut . Although we find 
that the data fit is almost independent of the value of 
T C ut we choose T cut = 0.005, which gives optimal coinci- 
dence with Bethe-ansatz results available for J2 = (see 
below). On the other hand, the ED data at very low tem- 
perature arc affected by finite-size effects. To circumvent 
this problem we proceed as follows. Wc first determine 
the temperature Ted down to which the first four dig- 
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FIG. 6: Specific heat calculated by RGM (solid) and ED 
(dashed curve; N=20) for J 2 = 0.15,0.18,0.2, and 0.24, from 
top to bottom. The inset exhibits the RGM results in an 
enlarged scale. Note that for J2 = 0.24 only ED data are 
shown. 
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FIG. 7: Finite-size dependence of the low-temperature spe- 
cific heat for J2 = 0.2 (upper panel) and 0.24 (lower panel). 
The lines represent ED data for N = 12, 14, 16, 18, 20 and 22, 
from top to bottom. 



its of the specific heat per site Cy(T) for N = 20 and 
N = 22 coincide. (We use the specific heat to determine 
Ted, because Cy(T) is most sensitive to finite-size ef- 
fects at low temperature, see also below.) Then we use 
the ED data points for xT 2 in the temperature region 
Ted < T < Ted + T cut to fit them to Eq. 0. We find 
that Ted varies from 0.22 at J 2 = to 0.03 at J 2 = 0.24. 
Obviously, we have to use ED data points at higher tem- 
perature for the fit in comparison to the RGM fit, in par- 
ticular at small values for J 2 . The results for y and y\ 
as well as for xg and X\ are shown in the insets of Figs. [3] 
and 2] It is obvious that the values for yo determined by 
RGM and ED are very close to each other. Note that for 
the unfrustrated ID ferromagnet the quantities yn and Xo 
were calculated by the RGM previously in Ref . [35|. It was 
found that y = 1/24 s=s 0.041667 and x = 1/4, which 
agrees with the Bethe-ansatz results of Refs. |32| and [33| 
[note that x defined in Ref. J32| is larger by a factor of 4 
than x given by Eq. ( Up], Our RGM data confirm these 
findings (see also Ref. [29|). The fitting of the ED data 
at J 2 = yields y = 0.0418, which is still in reasonable 
agreement with the Bethe-ansatz result. Including frus- 
tration J 2 > we find an almost linear decrease in j/o 
as well in xq with J 2 down to zero at J 2 = 0.25 (cf. the 
insets of Figs. [3] and [4}. A linear fit of the RGM data 
points yields the relations 

j/o = (1 - 4J 2 )/24 ; x = (1 - 4J 2 )/4 , (9) 

which describe the RGM data in high precision. The 
vanishing of z/o and of xq at J 2 = 0.25 reflects the 
zero-temperature phase transition at this point and in- 
dicates the change in the low-temperature behavior of 
the physical quantities at the critical point. Using the 
same J 2 data points as in the insets of Figs. [3] and fj] 
a polynomial fit according to y\ = a y + 6 y J 2 + c y J% 
(xi = a x + b x Ji + c x jf ), indeed, yields, at J 2 = 0.25, 
finite values y x = 0.047 for RGM and y x = 0.043 for ED, 
and x\ = 0.147 (RGM only). Hence, our data suggest 
a change in the low-temperature behavior of x (£) from 
X oc T~ 2 (£ oc T- 1 ) at J 2 < 0.25 to x oc T~ 3 / 2 (£ ex 
T -1 / 2 ) at the zero-temperature critical point J 2 = 0.25. 
Let us mention here again that our results for the critical 
indices 7 and v at J 2 = 0.25 are based on the validity of 
Eqs. (J7J) and (jSJ) and the extrapolation of our data from 
J 2 < 0.25 to J 2 = 0.25. A slightly different index 7 also 
being below the "ferromagnetic" value 7f = 2 discussed 
above, namely 7 = 4/3, is obtained^, if one employs the 
modified spin-wave theory by Takahashi^l at J 2 = 0.25. 
The next quantity we consider is the specific heat Cy . 
In Fig. [5] our RGM and ED results for Cy are compared 
with the TMRG data£ Obviously, the ED results are in 
a very good agreement with the TMRG data. The de- 
viation at low temperature, appearing for J 2 = 0.125 as 
an increased value of Cy for 0.02 < T < 0.1, is ascribed 
to finite-size effects (see also the discussion below). For 
larger values of J 2 the specific heat shows another in- 
teresting low-temperature feature (see Fig. [5]). In the 
region 0.125 < J 2 < 0.25 with a ferromagnetic GS, the 



specific heat exhibits two maxima. Besides the broad 
maximum at T s=s 0.6, an additional frustration-induced 
low-temperature maximum appears, which is found by 
the ED and RGM methods for J 2 > 0.125 and > 0.16, 
respectively. As shown by a detailed analysis (see also 
below), the behavior of Cy at very low temperature is 
appreciably affected by finite-size effects. In particular, 
in the ED data, the low-temperature maximum is super- 
imposed by a quite sharp extra finite-size peak, as can 
be clearly seen in Fig. [5] for J 2 = 0.24. In view of this, 
the height and the position of the true additional low- 
temperature maximum cannot be extracted unambigu- 
ously from the ED data, however, its existence is not 
questioned by this ambiguity. On the other hand, the 
RGM (see inset of Fig. yields a shift of the maximum 
to lower temperature with increasing frustration. 

To illustrate the finite-size effects at low temperature, 
in Fig.[7|the ED data for the specific heat for J 2 = 0.2 and 
0.24 and different chain lengths are plotted. As already 
discussed above, the first four digits of the Cy{T) data 
for N = 20 and 22 coincide down to T ED w 0.04 (T ED ~ 
0.03) for J 2 = 0.2 (J 2 = 0.24). (Note again that for J 2 = 
the corresponding value Ted ~ 0.22 is much larger.) 
Below Ted finite-size effects become relevant (cf. Fig.[7|). 
However, from Fig.[7]it is also evident that the extra low- 
temperature finite-size peak behaves monotonously with 
N. Hence a finite-size extrapolation of the height c pea k 
and the position T pea k of the extra peak is reasonable. 
We have tested several extrapolation schemes and found 
that a three-parameter fit based on the formula a(N) = 
«o + a\/N 2 + a,2 /N A is well appropriate to extrapolate 
both Cp ea k and T pea k to N — > 00. The results of such 
an extrapolation are shown as filled squares in Fig. [7] 
The extrapolated data points indicate that the extra peak 
indeed is a finite-size effect and it vanishes for N — > 00. 
However, it is also obvious that the characteristic steep 
decay of the specific heat down to T = starts at lower 
temperature T* when approaching the zero-temperature 
critical point (we find T* « 0.05,0.007, and 0.002 for 
J 2 = 0,0.2, and 0.24, respectively). This behavior is in 
accordance with the shift of the low-temperature RGM 
maximum in Cy mentioned above and is relevant for low- 
temperature experiments on quasi-lD ferromagnets. 

Finally, let us mention that in an early paper by Tonc- 
gawa and Harada£ and also recently by Heidrich-Meisner 
et al£ and Lu et al£ a double-maximum structure in Cy 
was already found for 0.25 < J 2 < 0.4, however, with 
a low-temperature maximum that becomes much more 
pronounced approaching the critical point. In this case, 
the low-temperature peak in Cy(T) was ascribed to ex- 
citations from a singlet GS to a low-lying ferromagnetic 
multiplet£ In our case J 2 < 0.25. Above the fully po- 
larized ferromagnetic GS multiplet many low-lying mul- 
tiplcts exist, and the appearance of the additional low- 
temperature maximum is attributed to a more subtle in- 
terplay between all of these low-lying states. 



V. SUMMARY 

In this paper we explored the influence of the NNN 
coupling J 2 < 0.25| Ji| on the thermodynamic properties 
of the ID spin-1/2 Heiscnberg ferromagnet using ED 
and RGM methods. The results of both methods 
are in qualitatively good agreement. We found that 
the critical behavior of the susceptibility x an d the 
correlation length £ is not changed by the frustrating 
J 2 . However, limT^o X? 2 and lim^^o £,T go to zero 
for J 2 — ► 0.25 1 Ji I indicating a change in the low- 
temperature behavior of x (£) from \ oc T -2 (£ oc T _1 ) 



at J 2 < 0.25| Ji| to x oc T" 3 / 2 (£ ex T" 1 / 2 ) at the critical 
point J 2 = 0.25|Ji|. Another interesting feature is the 
appearance of a double-maximum structure in the spe- 
cific heat Cy, where the additional frustration-induced 
low-temperature maximum was found by ED (RGM) to 
occur for J 2 /\ J x \ > 0.125 (0.16). 
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